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We compute numerically the effective potential for the (A$*)4 theory on the lattice. Three different methods were used to 
determine the critical bare mass for the chosen bare coupling value. Two different methods for obtaining the effective potential 
were used as a control on the results. We compare our numerical results with three theoretical descriptions. Our lattice data are 
in quite good agreement with the "Triviality and Spontaneous Symmetry Breaking" picture. 



1. Introduction. 

The response to an external field is a standard tool 
to obtain non-perturbative information in quantum 
field theories. This is even more important in a scalar 
field theory since the introduction of an external cur- 
rent gives directly the effective potential (up to an 
integration constant), thus providing important infor- 
mation on the stability of the system. 

The conventional interpretation of the triviality in 
(A<i>'*)4 theories is based on Renormalization-Group- 
Improved-Perturbation-Theory (RGIPT) |1-3|. This 
picture predicts a second-order phase transition and 
a vanishing Higgs mass nih in the continuum limit 
if V, the physical v.e.v, is held fixed. An alternative 
interpretation of triviality |@}]6l predicts a phase tran- 
sition that is very weakly first-order and that nih and 
V are both finite, cutoff-independent quantities. The 
latter interpretation originates |4l from the Gaussian 
approximation. The two alternative pictures can be 
distinguished by a sufficiently precise lattice calcula- 
tion of the effective potential V^ff to this end we per- 
form a model-independent "numerical experiment" to 
test the predictions of both the conventional RGIPT 
picture and the alternative picture. 

2. Effective potential on the lattice. 

To evaluate the effective potential on the lattice we 
introduce a constant external source J in the classical 
action for the (one-component) (A$'*)4 theory in Eu- 
clidean space. One of the possible lattice discretiza- 



tions reads 



^<i>2(a;) + ^<I'4(x)- J<i>(a;)| (1) 

where x stands for a generic lattice site, the lattice 
fields are expressed in lattice units (a — 1), and 
Ao > 0. For SSB the basic quantity is the expec- 
tation value of the bare scalar field ^{x) (B=Bare) 
($),/ = (f>B{J), since determining 4)b{J) at sev- 
eral J-values is equivalent to inverting the relation 
J — J{4>b) = dVeff/dipB involving the effective po- 
tential Veff{(pB)- In this way, starting from the ac- 
tion in Eq.([0, the effective potential of the theory is 
rigorously defined up to an arbitrary integration con- 
stant (usually chosen to fix Kff (0) = 0) and is convex 
downward |7 1. 

Since we want to simulate the (A$'*)4 lattice the- 
ory close to the continuum limit we have to deter- 
mine the critical value rc of the bare mass parameter 
j'o at J = 0, which defines the 'Coleman- Weinberg 
regime' where to, the mass gap of the symmetric 
phases, vanishes. The critical bare mass parameter 
Tc can be determined from the susceptibility x = 

1/f^iatt - Where $ = H^) 

is the average field for a given configuration and the 
brackets stand for the average on the lattice configu- 
rations produced in a Monte Carlo run. 

One expects that, near the critical region, ~ 
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{fc — ro), modulo logarithmic corrections to the free- 
field scaling law. One can thus determine rc by ex- 
trapolation to vanishing x~^- Strictly speaking, this 
method is valid only for a second-order phase tran- 
sition where the phase transition value tq = co- 
incides with Tc and such that, at rc — Vg, both m 
and rufi vanish. In the case of a very weak first- 
order phase transition where f4| \rc — rs\/rc ^ 
exp(— l/(2s)), s = 3Ao/167r^ ^ 1, the induced nu- 
merical uncertainty should be negligible. 

For the Monte Carlo simulation of the lattice field 
theory described by Eq. JO we followed the upgrade 
of the scalar field $(x) (using Metropolis) by the up- 
grade of the sign of $(a;) according to the embed- 
ded Ising dynamics (SJ. The zero external field spin- 
flip probability of the Swendsen-Wang algorithm is 
slightly modified to take into account the external cur- 
rent J. 

Our data fE\ are well described by the simple lin- 
ear fit — air — rd and do not show evidence 
of the logarithmic corrections. We evaluate the above 
quantity both in the symmetric and broken phase. We 
determined rc also through the following fit 1 9 1 to the 
generalized magnetization ($): 



($} ^ a{rc ~ ry^^\lTi\r ~ Tclf + 6 . 



(2) 



Combining the above estimates for rc we get rc = 
—0.2280(9), in perfect agreement with the indepen- 
dent analysis IIIOI which predicts, for Aq — 0.5 and 
L = 16, rc = -0.2279(10). Thus we have three 
independent and consistent evaluations of rc at Aq = 
0.5 on a IB"' lattice, that represent a precise input defi- 
nition of the 'Coleman- Weinberg regime' with the ac- 
tion Eq.lHJ at J 0. 

We have used two independent methods to compute 
the effective potential. Firstly, we ran simulations of 
the lattice action Eq. Q for 16 different values of the 
external source in the range 0.01 < | J| < 0.70. In 
this way, as outlined in Eqs. (2-4), we directly obtain 
the slope of the effective potential (from which Vcs 
can be obtained, up to an additive integration con- 
stant). We performed our numerical simulations by 
using for the sign upgrade both the Metropolis and 
the S-W cluster algorithm. 

As a third additional check of our results, we 
performed a calculation using an alternative ap- 
proach 111 II to V^ff based on the approximate effective 



potential t/cff ('/>_b; ^) defined through 

exp{-L/cff(0B;r2)} = 

J[D<^>] 6 (J^ Jd^x <^{x) - exp -5[$] . (3) 

In the limit in which the 4-volume ft ^ oo, UeS 
tends to the exact K;ff (0b). The difference between 
Ucs{4>B'-,^) and Vcs{4>b) gives both a consistency 
check of our calculations and an indication of the ef- 
fects due to the finiteness of our lattice. 

3. Comparing theory with the lattice data. 

We can now compare our three different sets of lat- 
tice data with the existing theoretical expectations. 

In the case of the picture in which triviality is com- 
patible with SSB |4 6| the predicted form (in the 
Coleman- Weinberg case, ro = rc, where no quadratic 
term is present in the effective potential) is: 







dVtriv 



a(t>% ln(0|) + 7(/)|, 



(4) 



where a and 7 are free parameters. (Their val- 
ues are approximation-dependent within the class of 
"triviality-compatible" approximations.) 

The RGIPT prediction exists in various slightly dif- 
ferent forms in the literature. We have first used the 
full two-loop calculation of Ford and Jones 1121 in 
the dimensional regularization scheme. The theoreti- 
cal prediction for J^^'™? (0b) = dt/2-loop/^^^ jjg. 
pends on two free parameters: the 't Hooft scale /i and 
the mass parameter of the classical potential. 

A different version of the RGIPT prediction, 
which re-sums various terms, is given in Eq. (242), 
Sect. 5.4.2, of the textbook |13| by Itzykson and 
Drouffe (ID), namely 



J 



ID/ 



In 



\'I>b\ 



1/3 



(47r)^0i 
18 In-' 



(5) 



and again we have two free parameters A and /i. 

The results of fitting the lattice data to the three the- 
oretical predictions, are shown in Table 1 . 

4. Conclusions. 

We have performed a numerical experiment to test 
two basically different and alternative pictures of 
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data 


jtriv 






Metropolis 


a = 0.0152(2) 

2 _ 15 
^ 16-2 


M = 8.0304(449) 

ivl — — U.UUZOl^lJ 

2 _ 142 
16-2 


l^il = 2.70(1) X 10^ 

/± — — U.UUOOl^oj 

2 _ 116 
16-2 


Swendsen-Wang 


a = 0.0152(1) 

2 _ 13 
^ 16-2 


M = 8.0128(283) 

M =— U.UUz4y(/J 

^2 _ 284 
16-2 


|/i| = 2.70(1) X 10^ 

/i = — U.UUoo(z ) 

^2 _ 223 
^ 16-2 


constraint eff. pot. 


a = 0.0156(2) 
7 = 0.4494(1) 

2 _ 10 
16-2 


fi = 7.9883(455) 
]VP = -0.0028(1) 

2 _ 109 
A 16-2 


1^1 = 2.69(1) X 10« 

A = -0.0063(3) 

2 _ 85 
A 16-2 



Table 1 

Results of the fits of J^™, j2-ioop^ ^j^^ jID ^^j. ^j^j-^g different sets of lattice data. 



'triviality' in the 'Coleman- Weinberg-regime' . Our 
results are striking and confirm that the lattice data 
cannot be reproduced by all theoretical models. In- 
deed the results of fitting the conventional theoretical 
predictions to our lattice data gives totally unaccept- 
able values of the x^/ f- On the other hand Eq. (0}, 
based on the alternative picture, gives excellent fits 
to all sets of data. This is an important evidence 
for the unconventional interpretation of "triviahty" of 
Refs. mil. 
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